Duty-ratio of cooperative molecular motors 
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Molecular motors are found throughout the cells of the human body, and have many different and 
important roles. These micro-machines move along filament tracks, and have the ability to convert 
chemical energy into mechanical work that powers cellular motility. Different types of motors are 
characterized by different duty-ratios, which is the fraction of time that a motor is attached to its 
filament. In the case of myosin 11 - a non-processive molecular machine with a low duty ratio - 
cooperativity between several motors is essential to induce motion along its actin filament track. In 
this work we use statistical mechanical tools to calculate the duty ratio of cooperative molecular 
motors. The model suggests that the effective duty ratio of non-processive motors that work in 
cooperation is lower than the duty ratio of the individual motors. The origin of this effect is the 
elastic tension that develops in the filament which is relieved when motors detach from the track. 

PACS numbers: 



I. INTRODUCTION 

Motor proteins are molecular machines that convert 
chemical energy into mechanical work by ATP hydrol- 
ysis. They "walk" on the microtubule and actin cy- 
toskeleton and pull vesicles and organelles across the cell 
Motor proteins can be classified into processive and 
non-processive motors. The former class includes motors 
like kinesins, which travel a long distance along their cy- 
toskeleton track (microtubules) without detachment [4. 
In the latter class, we find motors like myosins that make 
only a single step along their tracks (actin filaments) 
before disconnection [3|. While processive motors can 
move cargoes by operating individually, non-processive 
motors need to work in cooperation to generate substan- 
tial movement. Cooperative action of myosin motors is 
implicated in a variety of cellular processes, including the 
contraction of the contractile ring during cytokinesis [H , 
adaptation of mechanically activated transduction chan- 
nels in hair cells in the inner ear [5| , and muscle contrac- 
tion Q . Another important example of cooperative mo- 
tor dynamics is found in motility assays, where filaments 
lide over a surface densely covered by motor proteins 



One of the more interesting outcomes of cooperative 
action of molecular motors is their ability to induce bidi- 
rectional motion. Bidirectional motion is observed when 
a filament is subjected to the action of two groups of 
motors that engage in a "tug-of-war" contest and exert 
forces in opposite directions Q. The motor party that 
exerts the larger force determines the instantaneous di- 
rection of motion, which is reversed when the balance of 
forces shifts from one group to the other. Earlier theoret- 
ical models suggested that that the characteristic reversal 
time of the bidirectional dynamics (i.e., the typical dura- 
tion of the unidirectional intervals of motion) ,r r ev, grows 
exponentially with the number of motors N Q . In these 
earlier models, the moving filament was treated as a rigid 
rod which does not induce elastic coupling between the 



motors. A recent motility assay of myosin II motors and 
actin filaments with alternating polarities challenged this 
prediction. It was found that the characteristic reversal 
times of the bidirectional motion in this motility assay 
were macroscopically lar ge, but practically independent 
of the number of motors [lOJ . This observation has been 
explained by a model that accounts for the elasticity of 
the actin filament HH • It has been shown that the 
motors indirectly interact with each other via the ten- 
sile stress that they generate in the elastic filament. The 
elasticity-mediated crosstalk between the motors leads 
to a substantial increase in their unbinding rates, mak- 
ing each motor effectively less processive and eliminating 
the exponential growth of r rcv with N [l2j ■ 

The reduction in the duty ratio (the fraction of time 
that each motor spends attached to the filament track) 
of the cooperative motors can be explained as follows: 
During bidirectional motion, the elastic filament is sub- 
jected to a tug-of-war between motors that exert opposite 
forces, which leads to large stress fluctuations along the 
elastic filament. These stress fluctuations are the origin 
of the elasticity-mediated crosstalk effect. Detailed anal- 
ysis shows that the typical elastic energy stored in the 
actin filament scales as [10, [ill 
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where TV is the number of motors, n is the number of 
connected motors, and a is a dimensionless parameter 
(which is closely related to the parameter f3* , to be de- 
fined below in eq. EH- For actin- myosin II systems, a is 
rather small (for the motility assay described in ref. [l(| , 
a~2x 10~ 3 ); but the energy released upon the detach- 
ment of a single motor (n — > n — 1): AE/ksT = —aN, 
can be quite large if N is large. This last result implies 
that for large N, the transitions of motors between the at- 
tached and detached states are influenced by both ATP- 
driven (out-of-equilibrium) processes as well as by ther- 
mal (equilibrium) excitations. The latter will be domi- 
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FIG. 1: The actin elastic filament is represented as a chain of 
nodes, connected by identical springs with spring constant k s . 
Each node is either connected to a single myosin II motor - in 
which it experiences a force of magnitude fo, or disconnected 
- in which it experiences no force. 



Thus, h(x) gives the total force applied on the chain up 
to the point x, with h(x = 0) = and h(x = N) = n, 
where n is the number of monomers connected to mo- 
tors in a given configuration. The solid line in fig. [5] 
shows the function h(x) corresponding to the configura- 
tion of 5 nodes depicted in fig. [TJ To calculate the elas- 
tic energy of a configuration, we introduce the function 
gix) = h(x) — (n/N)x, which is depicted by the dashed 
line in fig. [2] and gives the total excess force accumulated 
up to x. The elastic energy can then be expressed as: 



nated by the changes in the elastic energy of the actin 
(see eq.[T]) and not by the energy of the individual motors. 

For the actin filament with alternating polarities, the 
elasticity-mediated crosstalk is a cooperative effect that 
reduces the degree of cooperativity between motors (de- 
creases r rov ) by decreasing their duty ratio (i.e., their 
attachment probability). But as discussed above, much 
of the strength of this effect is related to the large stress 
fluctuations that develop in the elastic filament due to 
the opposite forces applied by the antagonistic motors. 
In view of this fact, it is fair to question the significance 
of this effect in the more common situation where polar 
filaments move directionally under the action of a single 
family of motor proteins. In this work we analyze this 
problem and show that the elasticity-mediated crosstalk 
effect in this system is indeed much smaller, but not en- 
tirely negligible. Our analysis of this effect is presented 
in the following section |TTJ The magnitude of this effect 
in acto- myosin systems is estimated in section IIIII 



II. STATISTICAL-MECHANICAL MODEL 

In what follows we model the elastic actin filament as 
a chain of N equally spaced nodes connected by TV — I 
identical springs with spring constant k s (see fig. [T]) . In 
the chain reference frame, the i-th node is located at 
Xi = (i — I )Al, where i = 1, . . . , N and Al is the spacing 
between the nodes. For brevity we set Al = 1. The chain 
lies on a "bed" of motors, where each node may be either 
free and experience no pulling force (fi = 0), or attached 
to one motor in which case it is subjected to a force of 
magnitude fi — fo. Different configurations of the sys- 
tem are defined according to which nodes are connected 
to motors and which are not. For a given configuration 
j, the elastic energy of the chain is given by the sum of 
energies of the springs: = ^fj^ 1 F 2 /2k s , where Fi 
is the force applied on the i-th spring. The forces Fi are 
calculated as follows: We first calculate the mean force 
/ = (j2^ =1 fi) /N, and define the excess forces acting 

on the nodes: /* =/» — /• The force on the i-th spring 
is then obtained by summing the excess forces applied 
on all the monomers located on one side of the spring: 

p V** f * r* 

r i ~ Z^;=i J i — Z^;=i+i J i ■ 

It is more convenient to analyze the problem using con- 
tinuous functions. Let us introduce the function h(x) 
which, for Xi < x < ccj+i, has slope +1 if the monomer 
at Xi is connected to a motor, and a slope 0, otherwise. 
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where 
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is the ratio between the typical elastic energy of a spring 
/o/2fc s . and the thermal energy ksT. 

To determine the mean number of connected motors, 
one needs to calculate the partition function 
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where p is the attachment probability, i.e. duty ratio, 
of a single motor, and z n is the partition function of 
all the configurations with exactly n connected motors. 
The function z n can be calculated by tracing over all 
the functions g{x) corresponding to configurations with 
n connected motors. Mathematically, the condition that 
exactly n motors are connected can be expressed through 
the following constraint on the function h(x). 
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FIG. 2: The functions h(x) and g(x) which correspond to the 
configuration shown in fig. [T] 



3 



To allow an analytical solution, we approximate this constraint by setting a = 2, in which case eq. [5] can be 
expressed in terms of g{x) = h(x) — (n/N) x as 
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With eq. [6j the partition function z n is given by 



z n = B(n,N) fv[g(x)]exp 1-/3* f g 2 (x)dx] S 
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where S is Dirac's delta-function. The function B (n, N) is introduced in eq. [7] in order to compensate for the error 
introduced by the approximated constraint eq. |51 We will determine this function through the requirement that for 
13* = 0, i.e. in the absence of elastic crosstalk between the motors, 
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which is simply the number of ways to choose n out of N monomers. 

In order to calculate the partition function z n , we use the Fourier space representation of S(x), 



(8) 



S(x-a) = — I e w{x ~ a) dw, 
2m 



(9) 



and the Fourier series of g(x), 



N/2-1 
k=-N/2 



Substituting eqs. M an d EH into eq. [7] yields: 
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dw V [g k ] exp 
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Tracing over g k can be readily performed, giving 



x exp 



(10) 



z n — B (n, N) - — / dw exp 

2m ' 
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The integral over w; can be evaluated using the method of steepest descent, which yields: 

z n ~B{n,N) e G{w °\ 



(12) 
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where, 
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and wq satisfying 
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For p* < 1, one gets 
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From eqs. H [13 HH and [HI one finds that 
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Inserting eq. [T|)]into eq. 1141 and expanding G(wq) in powers of -\/]3~*, yields 



N/2 



(16) 



(17) 



G{ Wo )~G{w Q )\p.= -sJ n{N 2 n) /3*. (18) 
Finally, for /?* < 1, the partition function z n is obtained by substituting eqs. [171 and IT51 into eq. [JS] which gives: 
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In order to calculate the partition function Z, one needs to substitute eq. [15] into into eq. 2] which gives: 

N 



z = J2( N )p n ( 1 ~p) N - n ^p 



n(N - n) 
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In the thermodynamic limit (N 3> 1), the sum in eq. [201 is dominated by one term which corresponds to the mean 
number of attached motors (n) . This term is given by 



(n) = N 



p-(l-2p) A /^fi^ 



From eq. [5T] we identify the effective attachment probability as 



(21) 



(22) 



Notice that the second term on the right hand side of 
eq. [221 is anti-symmetric around p — 1/2, and that for 
p < 1/2 (p > 1/2), the effective attachment probability 
p e g is smaller (larger) than p. This observation is directly 
related to the fact the elasticity-mediated crosstalk effect 
is driven by the tendency to reduce the force fluctuations 
along the elastic filament. For p < 1/2 (p > 1/2) the 
force fluctuations are reduced by the detachment (attach- 
ment) of motors, which brings the system closer to the 
limiting case p = (p = 1) where the force fluctuations 
vanish. 

To test the validity and range of applicability of eq. [22] 
we conducted Monte Carlo (MC) simulations of elastic 
chains of N — 1000 monomers with p = 0.05, which 
is the typical duty ratio of myosin II motors [13j . Sys- 
tems corresponding to different values of /3* were sim- 
ulated using the parallel tempering method. The sim- 
ulations include two types of elementary moves (which 
are attempted with equal probability) - one in which 



the state (connected/disconnected) of a randomly cho- 
sen node changes, and the other in which two randomly 
chosen nodes with opposite states change their states si- 
multaneously. For each move attempt, the model en- 
ergy of the chain is recalculated, and the move is ac- 
cepted/rejected according to the conventional Metropo- 
lis criterion. Our MC results are summarized in fig. [3] 
For P* < 0.2, we find an excellent agreement between 
our computational results and eq. [55] (which has been 
derived for /3* <C 1). Notice that eq. [51] does not include 
any fitting parameters. For larger values of /3*, eq. 1221 
overestimates the decrease in p c g . 



III. ACTIN-MYOSIN II SYSTEMS 

Eq. [3] relates the dimensionless parameter f}* to three 
physical parameters of the system - the temperature T, 
the typical force exerted by the motors /o, and the effec- 
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FIG. 3: The effective attachment probability as a function 
of the dimensionless parameter f3* . The circles denote the 
results of the MC simulations, while the solid line depicts the 
analytical approximation for /3* <C 1, eg: 1221 

tive spring constant of the actin filaments segment be- 
tween between two binding sites of the motors, k s . The 
latter parameter can be further expressed as 

Y A 

K = — , (23) 

where Y is Young's modulus of the actin, A is the cross- 
sectional area of an actin filament and / is the distance 
between binding sites. For myosin II motors, forces in 
the range of fn = 5—10 pN have been measured ex- 
perimentally [HI • The actin-cross sectional area (in- 
cluding the tropomyosin wrapped around the actin helix) 
is A = 23 nm 2 , and the Young's modulus of the actin- 
tropomyosin filament is Y = 2.8 GPa [l5|,[l6j]. The value 
of I is somewhat more difficult to assess. One possibility 
is that I ~ 5.5 nm, which is simply the size of the G- 
actin monomers, each of which includes one binding site 
for myosin motors [ljj ■ Another possible value is related 
to the double helical structure of F-actin and the fact 
that it completes half a twist about every 7 monomers, 
i.e. every ~ 38 nm Since the binding sites follow 

a twisted spatial path along the double helix, many of 
them remain spatially unavailable to the motors. In the 
motility assay, the motors are located underneath the F- 
actin, and the distance between the binding sites along 
the line of contact with the bed of motors is I — 38 nm. A 
third choice for I, which may be more relevant to skeletal 
muscles, is I = 14 nm. This value is derived from the fact 
that within the sarcomere (the basic contractile unit of 
the muscle) an average number of three thick myosin fil- 
aments surround one thin actin filament. The separation 
between collinear motor heads along the thick filament is 
~ 43 nm [19r|2l| , and because the actin is surrounded by 
three thick filaments, the distance between the motors 
along the actin is I — 43/3 ~ 14 nm. 

For the duty ratio p, the range of experimental val- 
ues is scattered and varies from p = 0.01 — 0.02 [22j to 
p = 0.05 fl3j . The uncertainty may be partially related 



to the elasticity cross-talk effect discussed here. It also 
stems from the fact that p also depends on the distance 
and orientation of the motor head (both in space and 
time) with respect to the associated binding site. These 
may vary between the motors which, hence, should have 
different attachment probabilities |23j |. 

Using the above mentioned values of system parame- 
ters, one finds that for acto-myosin systems the parame- 
ter (3* lies within the range of 5 x 10~ 4 < (3* < 5 x 10~ 3 . 
Setting the duty ratio of myosin II to p = 0.05 (as used 
in the MC simulations) we find that for this range of j3*, 
the effective duty ratio is slightly lower than p and lies 
within the range of p c g = 0.97p (for /3* = 5 x 10~ 4 ) and 
Pcff = 0.90p (for /?* = 5 x 10~ 3 ). Using a lower estimate 
for the duty ratio p = 0.02 (2^|, we find that the effective 
duty ratio drops to 0.83p < p e g < 0.95p for the same 
range of f3* . 

As stated earlier, the cooperative action of myosin mo- 
tors "compensates" for the non-processive character of 
the individual motors. The force generated by a group of 
N motors is (F) = Np e gfo. Which force /o maximizes 
the effective force per motor f c g = (F)/N and, hence, 
the force production of the collectively working motors? 
From eq. [3] and [H] we find that the maximum value of 
/off is achieved when the force of the individual motors 
is 

/•max _ 2p k s kBT 

/o = T^-p)jW^pj (24) 

Setting the values of the system parameters as above 
(Y = 2.8 GPa, A = 23 nm 2 , I = 38 nm) and taking 
p = 0.02 as the duty ratio of a single motor, we find that 
jmax _ 25 pN; for which _ o.25 pN. For forces in 

the range of /o — 5—10 pN, which are typically mea- 
sured for myosin II motors [1 31 ] . the effective mean force 
per motor is f e g ~ 0.1 — 0.15 pN, which is about half of 
the optimal effective force f e g (/o nax ). We, thus, conclude 
that mysoin II motors work quite close to conditions that 
maximize their cooperative force generation. 

IV. DISCUSSION 

There has recently been a considerable interest in the 
collective behavior of molecular motors, especially in re- 
lation to cooperative dynamics of cytoskclctal filaments 
in motility assays. Most studies have focused on the bidi- 
rectional motion arising when the filament is driven by 
two groups of antagonistic motors, or in the case when 
one motor party works against an external force. The 
present work is motivated by our recent studies of bidi- 
rectional motion in acto-myosin motility assays which 
demonstrated that the duty ratio of the motors (and, 
hence, the level of cooperativity between them) is re- 
duced by the elasticity of the actin backbone. In this pa- 
per we extended our studies to motility assays in which 
similar motors act on a polar actin filament without a 
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counter external force. To single out the filament elas- 
ticity crosstalk from other possible collective effects (e.g., 
those associated with non-equilibrium ATP-assisted pro- 
cesses and with the elasticity of the motors themselves 
[24T - I261 ] h we neglect motor-to- motor variations and use a 
model in which the motors are characterized by two mean 
quantities: their attachment probability to a rigid (non- 
elastic) filament p, and the mean applied pulling force 
fo- We expect this mean field description of the motors 
to hold when the number of motors N becomes large. 
We calculate the attachment probability to the elastic 
filament p c g < p from the partition function associated 
with the filament elastic energy. The elastic energy of 
filament can be treated as an equilibrium degree of free- 
dom (of a system which is inherently out-of-equilibrium) 
because the mechanical response of the filament to the 
attachment/detachment of motors occurs on time scales 
which are far shorter than the typical attachment time of 
the motors (see discussion in ref. [T^ |. The assumption 
that the actin filament is in mechanical equilibrium is also 
made in theoretical studies of intracellular cargo trans- 
port H and muscle constraction (27j) Our calculation 
shows that p g is only slightly smaller than p. This re- 
sult is very different from our previous findings for motil- 
ity assays of bidirectional motion, where the elasticity- 
mediated crosstalk effect is substantial. Finally, we note 
that although in both this and previous studies we found 
that p e g < p [HI, HI] , the opposite relation cannot be ex- 
cluded under certain conditions (e.g., when the filament 
experiences external forces as in single molecule experi- 
ments or muscle contraction). We intend to address this 
opposite scenario in a future publication, in which we 
investigate the variations in both p (the ATP-hydrolysis 
related attachment probability) and p c g (which also in- 
cludes the contribution due to the elastic crosstalk effect) 
with the muscle contraction velocity. 

We thank Anne Bernheim and Sefi Givli for useful dis- 
cussions. 



motor forces /, = J™ otor — / i MF , we can rewrite eg. I All as 
fi - /r cous + F t - = 0. (A2) 



For the chain's center of mass (CoM), the equation of 
motion reads 
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0; (A3) 



but the viscous drag force is distributed uniformly along 
the chain (i.e., it is equal for all masses) and, therefore, 
eq. IA3I gives 



T N f 

rviscous Z-^z— 1 J 1 p 

Ji - N — j- 



Using this last result, eg. I A2 1 reads 
Fi - = -/*, 



(A4) 



(A5) 



where /* is the excess force. Since F = 0, we find that 
Fi = -A*. Then, F 2 = F 1 - f* = -f* - /£; and, in 
general, 



A' 
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Appendix A: The elastic energy 

We model the actin filament as a linear chain of iden- 
tical particles of mass m connected by elastic (massless) 
springs with spring constant k s (see fig. 0} . Three types 
of forces are exerted on the particles: (i) the motor forces 
/f lotor , (ii) the spring forces Fi, and (iii) friction drag 
forces /j dras . Because the motion is highly overdamped, 
the total instantaneous force on each mass vanishes. For 
the i-th mass, the equation of is 

yrmotar _ yrdrag + p. _ p.^ = Q ^ (M) 

with F = Fjy = 0. The drag force includes two contri- 
butions - one is due to motor friction (MF), / 4 MF , and the 
other due to friction with the viscous medium, jv iscous . 
The motor friction forces act only on the particles which 
are connected to the motors. Therefore, by redefining the 



spring forces — - - — — - - — — - - — 

ivi Ft Ft 

FIG. 4: The actin is modeled as a chain of masses and springs 
that moves in a highly viscous medium. Three types of forces 
are present in the system: the motor forces, a uniform viscous 
drag, and the spring forces. The motion is highly overdamped 
(zero acceleration) and, thus, the net force on each mass van- 
ishes. 

If the motion is only partially overdamped (including 
in the limit of no friction) , all the masses move together 
at the same acceleration. One can repeat the above cal- 
culation and show that eg. I A6I remains valid. 
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